Current and future ground-based interferometers require 
knowledge of the atmospheric time constant to, but this param- 
eter has diverse definitions. Moreover, adequate techniques for 
monitoring to still have to be implemented. We derive a new 
formula for the structure function of the fringe phase (piston) in 
a long-baseline interferometer, and review available techniques 
for measuring the atmospheric time constant and the shortcom- 
ings. It is shown that the standard adaptive-optics atmospheric 
time constant is sufficient for quantifying the piston coherence 
time, with only minor modifications. The residual error of a 
fast fringe tracker and the loss of fringe visibility in a finite 
exposure time are calculated in terms of the same parameter. A 
new method based on the fast variations of defocus is proposed. 
The formula for relating the defocus speed to the time constant 
is derived. Simulations of a 35-cm telescope demonstrate the 
, feasibility of this new technique for site testing. 
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Abstract. 



1. Introduction 

Astronomical sites for classical observations are character- 
ized in terms of atmospheric image quality (seeing). For high- 
angular resolution techniques such as adaptive optics (AO) and 
interferometry, we need to know additional parameters. The at- 
mospheric coherence time is one of these. Here we refine the 
definition of the interferometric coherence time, review avail- 
able techniques, and propose a new method for its measure- 
ments. 

The AO time constant, To, is a well-defined parameter re- 
lated to the vertical distribution of turbulence and wind speed 
(Roddier 1981). To correct wave fronts in real time, a suffi- 
cient number of photons from the guide star is needed within 
each coherence area during time to. This severely restricts the 
choice of natural guide stars and tends to impose the complex 
use of laser guide stars (Hardy 1998). It is shown below that 
new, simple methods of to monitoring are still needed. 

Modern ground-based stellar interferometers attain ex- 
treme resolution, but their sensitivity is limited by the atmo- 
sphere. Even at the best observing sites, such as Paranal in 
Chile, fast fringe tracking is not fully operative yet, and one 
therefore tends to employ exposure times that are short enough 
to "freeze" the atmospheric turbulence. The price is a substan- 
tial loss in limiting magnitude. It is hence important to measure 
the time constant, to, of the piston - i. e. the mean phase over 
the telescope aperture - at existing and future sites. However, 
the exact definition of to is not clear, any more than are methods 
to measure it. Do we need an interferometer to evaluate fo? Is to 
different from to? Does it depend on the aperture size and base- 
line? We review various definitions of the interferometric time 
constant based on the piston structure function (SF), on the er- 
ror of a fringe tracker, and on the loss of fringe contrast during 
a finite exposure time. It is shown that the piston time constant 



is proportional to the AO coherence time To, both depending on 
the same combination of atmospheric parameters. 

During site exploration campaigns, one would like to pre- 
dict the performance of large base-line interferometers, and it is 
desirable to do this with single-dish and, preferably, small tele- 
scopes. The existing techniques for to measurement are listed 
and a new method for site testing proposed. 

2. Atmospheric coherence time in interferometry 

2.1. Atmospheric coherence time t 

First, we introduce the relevant atmospheric parameters and the 
AO time constant to. For convenience, we outline the essential 
formulae, but for the general background, we refer the reader 
to Roddier SlW\\ . 

The spatial and temporal fluctuations of atmospheric phase 
distortion ip are usually described by the SF 

Zyr, f) = ([cp(r', f) - tfr + r',t+ t')f), (1) 

which depends on the transverse spatial coordinate r and time 
interval t. The angular brackets indicate statistical average. 

The atmosphere consists of many layers. The contribution 
of a layer i of thickness dh at altitude h to the turbulence in- 
tensity is specified in terms of Cl(h)Ah, equivalently expressed 
through the Fried parameter r^ 5 J = 0.423k 2 Cl(h)dh, k — 2k /A 
being the wavenumber. The spatial SF in the inertial range (be- 
tween inner and outer scales) is 

£V(r,0) = 6.883 (|r|/r ) 5/3 . (2) 

It is assumed that each layer moves as a whole with the 
velocity vector V(/i) (Taylor hypothesis). The temporal SF of 
the piston fluctuations D^j(0, t) in one small aperture due to a 
single layer is then equal to, the spatial SF at shift Vt, 
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ZV/(0, f) = 6.883 [V(h)t/r ,i] 



(3) 
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Summing the contributions of all layers, we obtain 

D v (0,t) = 2.910 f 5/3 k 2 V 5/3 (h)C 2 n (h)dh 
Jo 

= 6.883 (fy 5 /3M,) 5/3 = (f/T ) 5/ \ 



(4) 



where To = 0.314 (V0/V5/3) is the AO time constant (Roddier 
1981) and the average wind speed V p is computed as 



V P = 



VP(h)C 2 „(h)dh 



Up 



(5) 



The formulae are valid for observations at zenith. At angle 
y from the zenith, the optical path is increased in proportion 
to the air mass, sec y, and the SF increases by the same fac- 
tor. Further, the transverse component of the wind velocity 
changes. In the following, we neglect these complications and 
consider only observations at zenith, but the analysis of real 
data must account for y + 0. 

2.2. Piston time constant 

In an interferometer with a large baseline (B » Lq, where Lq: 
turbulence outer scale) the phase patterns over the apertures are 
uncorrelated on short time scales. Thus, for a small time inter- 
val (f < B/V), the SF of the phase difference <p (do not confuse 
with the phase if) in an interferometer with two small apertures 
will simply be two times larger, D^it) - 2D^(0, t) (Conan et 
al. I1995i . As a result the differential piston variance reaches 1 
rad 2 for a time delay to = 2~ 3 ^ 5 tq = 0.66 to. Note that, in the 
case of smaller baselines and large outer scales - when the as- 
sumption B » Lq becomes invalid - D^t) < 2D^(0, t) and the 
resulting coherence time, accordingly, lies between 0.66 To and 
To- Yet, B » Lq applies to the characterization of large baseline 
interferometers at low-turbulence sites. 

When an interferometer with larger circular apertures of 
diameter d is considered, phase fluctuations are averaged in- 
side each aperture. As shown later, for time increments smaller 
than d/V, the piston structure function is quadratic in t and is 
essentially determined by the average wave-front tilt over the 
aperture. The variance of the gradient tilt a (in radians) in one 
direction is (Roddier 198 lj Conan et al. [l995l Sasiela fT994l l 



0.170 A 2 r Q 5/3 d- 1/3 . 



(6) 



We write the piston SF in this regime as D,p{t) as 2 {kcr a Vf) 2 , 
sum the contributions of all layers, and obtain the expression 



£V(?) as 13.42 (V 2 t/r ) 2 (r /d) 



1/3 



(tlhf 



(7) 



>l/6 



where the modified time constant t\ = 0.273 (ro/Va) (d/ro) 
The analysis of the tilt variance with finite outer scale by Conan 
et al. (2000 1 is applicable here. The finite outer scale reduces 
the amplitude of the tilt and hence increases the piston time 
constant, but this effect depends on the aperture size and is not 
very strong for d < 1 m. 

Note that for small time intervals there is a weak depen- 
dence of the SF on the aperture diameter. Also, the wind ve- 
locity averaging is slightly modified. However, the expressions 
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Fig. 1. Theoretical temporal power spectrum of the fringe posi- 
tion at 0.5 fim wavelength. The two telescopes are separated by 
100 m and have mirrors of 2 m diameter, the Fried parameter 
equals ro=ll cm, the wind vector makes an angle of a — 45° 
with the baseline, V - lOm/s. The vertical lines correspond 
to the frequencies: 0.2 V/B and 0.3 V/d. The asymptotic power 
laws are v~ 2 / 3 , v~ 8 / 3 , y~ 17 / 3 from lowest to highest frequencies. 



for t\ and to produce similar numerical results as long as d/ro 
is not too large. Thus, the system-independent definition of the 
AO time constant ® also gives a good description of the tem- 
poral variations of the piston. 

For time delays of approximately B/V and larger, the pis- 
tons on two apertures are no longer independent. However, es- 
timates of the time interval over which the Taylor hypothesis 
is valid range from ~ 40 ms (Schoeck & Spillar [T998> to sev- 
eral seconds (Colavita et al. 1987 1. Hence, at time intervals of 
1 s or more, the Taylor hypothesis is insecure. Moreover, the 
finite turbulence outer scale reduces the amplitude of slow pis- 
ton variations substantially. Here we concentrate only on rapid 
piston variations where our approximations are valid. 

2.3. Piston power spectrum and structure function 

The temporal power spectrum of the atmospheric fringe posi- 



tion has been derived by Conan et al. J1995l >. Their result is 
reproduced in Appendix A with minor changes. The temporal 
piston power spectrum (IA.4> produced by a single turbulent 
layer is represented in Fig^for a specific set of parameters. 
Because of the infinite outer scale Lq, this example is not re- 
alistic for frequencies below ~ 1 Hz. Moreover, as discussed 
in Sect. 12.21 Taylor's frozen flow hypothesis becomes invalid 
at low frequencies. Due to the infinite Lq, the asymptotic be- 
havior of the spectrum, and in particular the cut-off frequen- 
cies, do not depend on the wind direction (Conan et al. 1995 1, 
whereas, in the real case of a finite outer scale, the cut-off fre- 
quencies are affected by wind direction, as described by Avila 
et al. d!997l . Conan et al. ( 1995 1 point out that changing turbu- 
lence intensity and wind speed shift the spectrum vertically and 
horizontally, respectively, without changing the shape of the 
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Fig. 2. Relation between average wind velocities V5/3 and V2 
for 26 balloon profiles at Cerro Pachon in Chile (Avila et al. 
2000 1. The full line corresponds to equality, the dashed line is 
V 2 = l.lVs/3. 



40 F 




t, ms 

Fig. 3. Structure function of the fringe position for an inter- 
ferometer with mirror diameters d = 0.1m, ro = 11cm, 
V = lOm/s. The vertical line corresponds to t — d/V. For 
t < d/V, the SF is quadratic in t (dotted line), cf. Eq.0 For 
longer time scales, Dp ~ 2D V (dashed line). 

curve on the log-log plot. In observations with a small baseline 
(~ 12 m), the proportionality to y~ 2 ^ 3 at low frequencies and to 
v~ 8 ^ 3 at medium frequencies has actually been measured, e.g. 
by Colavita et al. JT91T71 . 

Based on the piston power spectrum, we derive in 
Appendix A the new expression of the piston SF valid for time 
increments t < min(B/V,Lo/V): 

D,(0 * 13.76 (V?M)) 2 [1.17 (d/r ) 2 + (Vt/r ) 2 ]- 1/6 . (8) 

As seen in Fig.[3] for t > d/V, the piston averaging over aper- 
tures is not important and we obtain = 2D V in agreement 
with heuristic arguments. For very short increments t <sc d/V, 
reduces to |7}. The average wind speed is V ~ V5/3 ~ V2. 
The difference between V5/3 and V2 is indeed small (Fig.|2}. 
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Fig. 4. Variance of corrected fringe position as a function of the 
bandwidth frequency of the correction system. The parameters 
of the simulation are identical to those of Fig.^ At frequencies 
higher than v c - 0.3V/d (vertical line), the variance is approxi- 
mated by (2nv c t\Y 2 (dotted line). 

2.4. Error of a fringe tracking servo 

A fringe tracker measures the position of the central fringe and 
computes a correction. The actual compensation equals the in- 
tegrated corrections applied after each iteration. Our analysis is 
similar to the classical work by Greenwood & Fried ( 1976 ). For 
a more detailed model that takes the effect of the finite expo- 
sure and response times of the phasing device into account, see 
the work by Conan et al. ( 2000b 1. The error transfer function of 
a first-order phase-tracking loop equals 

T(v) = iv/(v c + iv), (9) 

where v c is the 3 dB bandwidth of the system. The temporal 
power spectrum of the corrected fringe position is w c (v) = 
|r(v)| 2 w^(v). The residual piston variance characterizes the per- 
formance of the phasing device. This variance is shown in 
Fig-E]as a function of v c and is given by 

v 2 /(v 2 + v 2 ) vi>(v)dv. (10) 

00 

When v c < Q3V/d, the fringe tracker is too slow and leaves 
a large residual error; only fast trackers with v c > Q.3V/d 
are of any practical interest. In this case, the dominant con- 
tribution to the residual variance in fllOi comes from the fre- 
quencies just below 0.3V/ d, where the filter is approximated 
as (v/v c ) 2 . Hence the residual variance is proportional to the 
variance of the piston velocity. There is a simple relation be- 
tween the residual error of the fringe tracker and the structure 
function of the piston. For small arguments t, we can replace 
2[1 - cos(27rv?)] * (2nvt) 2 in the expression JA.5t for the phase 
SF. Then the residual error of the fast fringe tracker is simply 

(r 2 (v c ) * LV[1/(2ttv c )] * {2nv c hY 2 . (11) 

Thus, we have established that the error of the fast fringe 
tracker and the initial quadratic part of the piston SF are es- 
sentially determined by the variance of piston velocity which, 
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in turn, depends on the tilt variance and the average wind speed 

v 2 . 

2.5. Summary of definitions and discussion 

Tabled assembles different definitions of the atmospheric co- 
herence time. We have demonstrated that the time constant to 
of the piston SF is proportional to the AO time constant To. For 
small time increments, a slightly modified parameter t\ should 
be used. 

A different, but essentially equivalent, definition of the pis- 
ton coherence time Tq = 0.81 ro/V5/3 = 2.58 to has been 
given by Tango & Twiss ( 1980) and reproduced by Colavita 
et al. J1987l >. It is the integration time during which the piston 
variance equals 1 rad 2 . When fringes are integrated over a time 
To, the mean decrease in squared visibility equals l/e. Here 
we use the more convenient definition ?o = 0.66 To based on 
the temporal SF and warn against confusion with Tango's Tq. 
The definition of Tq is valid only for T > d/V, while shorter 
integration times are of practical interest (see below). 

The performance of the fringe-tracker in a long-baseline 
interferometer can be characterized by the atmospheric time 
constant t\ or, equivalently, by the average wind speed Vi- The 
AO time constant to (or V5/3) is also a good estimator of the 
piston coherence time, especially for small apertures d ~ ro. 

In order to reach a good magnitude limit, all modern inter- 
ferometers have large apertures d > tq. The atmospheric vari- 
ance over the aperture is 1.03 (t//r ) 5 ' 3 > 1 rad 2 and has to be 
corrected by some means (tip-tilt guiding, full AO correction, 
spatial filtering of the PSF) even at short integration times. The 
temporal piston variance will also be >1 rad 2 on time scales of 
approximately ro/V and longer. Hence exposure times shorter 
than ro/V or fast fringe trackers are required in order to main- 
tain high fringe contrast. In this regime, the relevant time con- 
stant that determines the visibility loss is t\, rather than To and 
T . 

All definitions of atmospheric time constants contain a 
combination of ro and V. As turbulence becomes stronger, the 
time constant decreases, although the wind speed may remain 
unchanged. Being less correlated, the parameters ro, V are thus 
more suitable for characterizing atmospheric turbulence than 
the parameters ro, To. Astronomical sites with "slow" or "fast" 
seeing should be ranked in terms of V rather than to. A fair 
correlation between V and the wind speed at 200 mB altitude 
has been noted by Sarazin & Tokovinin i2002t . 

3. Measuring the atmospheric time constant 

3.1. Existing methods of t measurement 

Table|2]lists methods available for measuring the atmospheric 
coherence time to or related parameters. The 3rd column gives 
an indicative diameter of the telescope aperture required for 
each method. Short comments on each technique are given be- 
low. 

SCIDAR (Scintillation Detection And Ranging) has pro- 
vided good results on to. It is not suitable for monitoring be- 
cause manual data processing is still needed to extract V(h), 
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Fig. 5. Five consecutive ring images distorted by turbulence 
and detector noise. Each image is 16x16 pixels (13.8"), the av- 
erage ring radius is 3", the interval between images is 3 ms, the 
wind speed is lOm/s. 

despite efforts to automate the process. Balloons provide only 
single-shot profiles of low individual statistical significance. 
The AO systems and interferometers give reliable results, but 
are not suitable for testing new sites or for long-term monitor- 
ing. 

The methods listed in the next four rows of Table[2]all re- 
quire small telescopes and can thus be used for site-testing. 
However, all these techniques have some intrinsic problems. 
SSS (Single Star SCIDAR) essentially extends the SCIDAR 
technique to small telescopes: profiles of C 2 (/z) and V(h) are 
obtained with lower height resolution than with the SCIDAR, 
and are then used to derive the coherence time. The GSM 
(Generalized Seeing Monitor) can only measure velocities of 
prominent layers after careful data processing. A coherence 
time, taa - which, however, does not have a similar depen- 
dence on the turbulence profile than to and r 1 — is deduced 
from the angle of arrival fluctuations. MASS (Multi- Aperture 
Scintillation Sensor) is a recent, but already well-proven, turbu- 
lence monitor. One of its observables related to scintillation in a 
2 cm aperture approximates V5/3 (Tokovinin 2002), but this av- 
eraging does not include low layers and thus gives a biased esti- 
mate of to. An even less secure evaluation of to can be obtained 
from DIMM (Differential Image Motion Monitor) by combin- 
ing the measured ro with meteorological data on the wind speed 
(Sarazin & Tokovinin l2002> . 

We conclude from this brief survey that a correct yet simple 
technique for measuring To with a small-aperture telescope is 
still lacking. Such a method is proposed in the next section. 

3.2. The new method: FADE 

To measure the interferometric or AO time constant, we need 
an observable related to V2 or V5/3. The atmosphere consists 
of many layers with different wind speeds and directions, so a 
true C 2 -weighted estimator is required. Its response should 
be independent of the wind direction. 

Wavefront distortions are commonly decomposed into 
Zernike modes (Noll 1976). The first mode, piston, cannot be 
sensed with a single telescope and the two subsequent modes, 
tip and tilt, tend to be corrupted by telescope vibrations. Of the 
remaining modes, the next three - defocus and two astigma- 
tisms - have the highest variance and are the best candidates 
for measuring atmospheric parameters. 

The total turbulence integral (or ro) is typically measured 
by the DIMM (Sarazin, & Roddie r[T990l . Lopez ( fT^92l tried to 
derive To from the speed of the DIMM signal, but this method 
did not prove to be practical. Because of its intrinsic asymme- 
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Table 1. Definitions of atmospheric time constants 



Quantity of interest Formula Time constant 



Phase SF 


D p (i) = 


a/T ) 5/3 


To 


= 0.314 ro/Vs/3 


Piston SF, t < d/V 


Df(t) = 


(t/hf 


h 


= 0.273 (r /V 2 )(d/r Q ) [ / 6 


Piston SF, t > d/V 


D^t) = 


(t/k) 513 


to 


= 0.66 r 


Piston variance during an exposure t > d/V 


^(0 = 


(t/T ) 5 ' 3 


To 


= 2.58 t 


Phase tracker error, v c > 0.3 d/V 


0>c) = 


= {27!V c h)- 2 


h 





Table 2. Methods of tq measurement 



Method 


Measurables 


d, m 


Problems 


Reference 


SCIDAR 


C 2 „(h), V(h) 


>1 


Needs large telescope 


Fuchs et al.rEMl 


Balloons 


ClQi), V(h) 


none 


Expensive, no monitoring 


Azouit & Verninl2005l 


AO system 


ro, t 


>1 


Needs working AO 


Fusco et al.l2004l 


SSS 


ClQi), V(h) 


>0.4 


Low height resolution 


Habib et al.l2006l 


GSM 


ro, V, T AA 


4x0.1 


No obvious relation to ro and t\ 


Ziad et all2000l 


MASS 


T o 


0.02 


Biased (low layers ignored) 


Kornilov et al.l2003l 


DIMM 


'"() 


0.25 


Indirect t estimate 


Sarazin & Tokovinin 2002 


FADE 


ro, t\ 


0.35 


New method 


This work 




try, DIMM does not provide an estimator of V that is inde- 
pendent of the wind direction. On the other hand, the fourth 
Zernike mode (defocus) is rotationally symmetric. 

We show in Appendix B that the variance of defocus ve- 
locity provides an estimator of the time constant t\ . The vari- 
ance of the defocus itself gives a measure of ro. Thus, we can 
measure both ro and The method is based on series of fast- 
defocus measurements, and we call it FADE (FAst DEfocus). 
The details of the future FADE instrument still need to be 
worked out and will be a subject of the forthcoming paper. Here 
we present numerical simulations to show the feasibility of this 
approach. We simulated a telescope of d — 0.35 m diameter 
with a small central obstruction e = 0.1. A conic aberration 
was introduced to form ring-like images (Fig.|5j. This configu- 
ration resembles a DIMM with a continuous annular aperture. 
The ring radius 3" was chosen. 



Monochromatic (A = 500 nm) images were computed on a 
64 2 pixel grid from the interpolated distortions and binned into 
CCD pixels of 0.86" size. We simulated photon noise corre- 
sponding to a star of R = 2 magnitude and 3 ms exposure time 
(20000 photons per frame) and added a readout noise of 15 
electrons rms in each pixel. 

The radius p of the ring image is calculated in the same 
way as standard centroids, by simply replacing coordinate with 
radius. The radius fluctuations Ap serve as an estimator for the 
defocus coefficient a\. The radius change is approximated by 
the average slope of the Zernike defocus between inner and 
outer borders of the aperture: 

Ap = C p a 4 ~ [2 V3(l + e)/n (A/d)] a A . (12) 

The complex amplitude of the light distorted by two phase 
screens at and 10km altitude with combined r ( > = 0.1 m was 
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pre-calculated on a large square grid (15 m size, 0.015 m pix- 
els). This distribution is periodic in both coordinates, and it was 
"moved" in front of the aperture in a helical pattern with the 
wind speed V to simulate the temporal evolution of the wave- 
front. The exposure time At = 3 ms corresponds to a wave- 
front shift VAt = 0.06 m for V = 20m/s, such that the initial 
quadratic part of the defocus SF (J3 = 2Vt/d < 1) extends only 
to ~ 3 At. 

Figure|6]shows the structure function, D p , of the ring-image 
radius calculated from several seconds of simulated data. It 
contains a small additive component due to the measurement 
noise (in this case 0.05" rms), which was determined from the 
data itself by a quadratic fit to the 2nd and 3rd points and its 
extrapolation to zero. The dashed lines are the theoretical SFs 
of defocus computed by (IB.5> and converted into radius with 
the coefficient C p H2\ . The slope between the second and third 
points of the simulated SF closely matches the analytical for- 
mula. 

To measure the speed of defocus variations, it is sufficient 
to fit a quadratic approximation to the initial part of the mea- 
sured SF, D p (t) « at 2 . Considering the noise, the best estimate 
of the coefficient a is obtained from the second and third points, 
a = [D p (2At) - D p (At)]/(3At 2 ). This estimator is not biased by 
white measurement noise. Equating the quadratic fit to the the- 
oretical expression D p (t) = 0.0269 (C p t/ti) 2 , we get a recipe 
for calculating the time constant from the experimental data, 

f, ~ 0.284 C p At [D p (2At) - D p (At)]- 1/2 . (13) 

Application of this formula to the simulated data gives 1 1 values 
of 3.88 and 2.20ms for wind speeds 10 and 20m/s, while the 
input values are 3.36 and 1.68 ms. Our simulated instrument 
slightly over-estimates t\ because the chosen exposure time of 
3 ms is too long. Indeed, the error gets worse for a higher wind 
speed and disappears for V — 5 m/s (true and measured t\ are 
6.73 and 6.62 ms) or for a shorter exposure time. In the real sit- 
uation of a multi-layer atmosphere, the experimental SF will be 
the sum of the SFs produced by different layers. The contribu- 
tion to the "jump" of the SF D p (2At) - D p (At) from fast layers 
will be reduced (in comparison with the quadratic formula) and 
will cause a bias in the measured t\, increasing its value. 

The crudeness of our simulations (discrete shifts of the 
phase screen, approximate C p , etc.) also contributes to the mis- 
match. Averaging of the image during finite exposure time has 
not been simulated yet. The response and bias of a real instru- 
ment will be studied thoroughly by a more detailed simulation. 
However, the feasibility of the proposed technique for measur- 
ing t\ is already clear. 

The next two Zernike modes number 5 and 6 (astigmatism) 
are not rotationally symmetric. However, the sum of the vari- 
ances of the velocities of two astigmatism coefficients is again 
symmetric. In fact, it has the same spatial and temporal spectra 
as defocus, with a twice larger variance. Therefore, simulta- 
neous measurement of the two astigmatism coefficients can be 
used to estimate the atmospheric time constant in the same way 
as defocus. Other measurables that are symmetric and have a 
cutoff at high frequencies can be used as well. However, defo- 
cus and astigmatism have the largest and slowest atmospheric 
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variances making it easier to measure than other higher-order 
modes. 

The FADE technique can be applied in a straightforward 
way to the analysis of the AO loop data, as a simple alterna- 
tive to the more complicated method developed by Fusco et al. 

num . 

4. Conclusions 

We reviewed the theory of fast temporal variations in the phase 
difference in a large-baseline interferometer. For a practically 
interesting case of large apertures d > ro, the piston SF usually 
exceeds 1 rad 2 at the aperture crossing time t = d/V. Hence, 
shorter times are of interest where the piston SF is quadratic 
(rather than oc t 5 ^). The relevant atmospheric time constant is 
1 1 . However, the standard AO time constant ro also provides a 
good estimation of the piston coherence time. Both these pa- 
rameters essentially depend on the turbulence-weighted aver- 
age wind speed V. 

A brief review of available methods for measuring ro shows 
the need for a simple technique suitable for site testing or mon- 
itoring, i.e. working on a small-aperture telescope. The FAst 
DEfocus (FADE) method proposed here fulfills this need. We 
argue that, for a given aperture size, this is the best way of ex- 
tracting the information on To. The feasibility of the method is 
proven by simulation, which opens a way to the development 
of a real instrument. An instrument concept using a small tele- 
scope, some simple optics, and a fast camera will be described 
in a subsequent article. 

Appendix A: Derivation of the piston structure 
function 

The spatial power spectrum of the piston is derived from the 
spatial atmospheric phase spectrum (Roddier ll981> 

(f +L 2 r I1/6 C\ dh, (A.l) 

o 

where f is the spatial frequency, Lq the turbulence outer scale at 
height h, and the other notations were introduced in Sect. 12.11 
We drop the explicit dependence of C„, Lq, and all following 
altitude dependent-parameters on h, to ease the reading of the 
formulae. The spatial filter that converts W v (£) into the piston 
power spectrum W^(f) is 

M 2 (f) = [2 sin(TrfB) A(f)] 2 (A.2) 
W^f) = M 2 (f) W„(f), (A.3) 

for a baseline vector B and the aperture filter function A(f). For 
a circular aperture of diameter d, A(f) = 2J\{nfd)l{nfd) and 
/ = |f |. There J„ stands for the Bessel function of order n. 

As usual, we assume that turbulent layers are transported 
with wind speed V directed at an angle a with respect to the 
baseline. The temporal power spectrum of the piston is then 
obtained by integrating in the frequency plane over a line dis- 
placed by f x = v/V from the coordinate origin and inclined at 
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angle a. Let f y be the integration variable along this line and (cf. Eq. 20 in Noll l 19761 . It follows that for t > d/V 
f 2 -fx + fy - The temporal spectrum equals 



D (f) * 13.77 (V5/3 tin,) 513 = (f/f ) 5/3 . 

v I W \f x cos a + f y sin a, f. cos or - f x sin or) d/ y For t < d/V, K\ (B) * 0.864 y3 2 and 

, f + °° , , f + °° , , ,,,, DM) ~ 13.41 (V 2 t/r ) 2 (r /d) 1/3 = (t/h) 2 

0.0388 £ 2 y-'C 2 d/! (/ 2 + L n 2 r I1/6 

J() J-oo 



(A.ll) 



(A. 12) 



We recover (7). This proves that the initial part of the piston SF 
[sin (nBf x cos a + nBf y sin a) A(/)] Z d/ v . (A.4) is indeed defined by the overall wavefront tilts. 

For a single turbulent layer, the piston SF is directly propor- 
tional to K\(B). Considering the small difference between two 
alternative definitions of the average wind speed, V5/3 as V 2 ~ 
V, a good approximation for the SF at all time increments will 
be 



We use the rotational symmetry of the aperture filter. This for- 
mula can be found in Conan et al. J1995I in a slightly dif- 
ferent form. The temporal power spectrum is defined here on 
v = (-00, +00) to keep the analogy with spatial power spectra. 
The temporal structure function of the piston is 

/-»+oo 

DJt) = 2[1 -cos(27rfv)] wJv) dv. (A.5) 



For an interferometer with a large baseline B » d, the width of 
the aperture filter is much larger than the period of the sin 2 fac- 
tor in iAA\ . We can then replace the sin 2 with its average value 
0.5. Assuming also that Lq » d, we obtain an approximation 
for the piston power spectrum 



w (v) * 0.0194 k 2 



OO 



V" 1 C 2 „ dh 



A z (/)/ -ll/3 dfy 



(A.6) 



With this approximation, 

D^t) = 0.0388 k 2 I C 2 d/i I I [1 - cos(2mf x V)] 

Jo J J-00 

x A 2 (/)r 11/3 df x df y 

C\dh I [1 -_/ (27rfV/)] 
Jo 



x A 2 (f)f- 8 ' 3 df 



(A.7) 



We used the relation (Gradshteyn & Ryzhik 1965 1 : 
J cos(27rz cos 8)d0 = 2n J${2nz). For a circular aperture of 
diameter d, 



D4.it) 



1.641 kf-d 



2^5/3 



f 

Jo 



C 2 d/; K x (2tV/d), 



(A.8) 



where the new dimensionless variables are B = 2tV/d and x 
nfd and the function K\ (B) 



Ki(J3) 



f 

Jo 



1.1183 



[2Ji(x)/x] V 8/3 [1 - 7 (y8x)] dx 
/? 2 



(A.9) 



(4.7 +B 2 ) 1 / 6 ' 

The approximation of K\{B) is accurate to 1% for all values 
of the argument and reproduces the analytic solutions of the 
integral for very large and very small B. For example, for large 
B the aperture filter tends to one; hence 



K X {B~) * x~ 8/3 [1 - /oOfo)] dx 
Jo 



= 7r/[2 8/3 r 2 (ll/6)sin(57r/6)]/? /j = 1.1 183^(^10) 



5/3 



9V3/ 



(A.13) 



D^t) * 3.88 (d/r ) 5/3 K x (2tV/d). 
With the approximation jA.9> . we finally obtain (JSJ. 

Appendix B: Fast focus variation 

The temporal power spectrum of the Zernike defocus coeffi- 
cient «4 is given in Conan et al. {1995 1 as 



w 4 (v) = 0.00969 ^ 



f 



CO 



y _1 c 2 d/z 



A^(/).r 11/3 d/ v 



(B.l) 



where A^(f) = 2 y/3Ji(nfd) / (nfd) is the spatial filter corre- 
sponding to the defocus on a clear aperture of diameter d (Noll 
[E>76l l. = v/V, f 2 = fl + / 2 , and we assume U » d. This 
expression is similar to JA.6> but has a two times smaller coef- 
ficient and a different aperture filter. The variance of defocus is 
a function of the Fried parameter: 



01 



X+CO 
W4(v)dv 
CO 

y-»+CO S\ y-»+CO 

= 0.00969 k 2 I C 2 d/z I I A 

Jo J J-00 



4(/)/ 



— 11/3 



= 0.0232 (d/r ) 



5/3 



df x df y 
(B.2) 



The variance of the defocus velocity has the following depen- 
dence on atmospheric parameters: 



CO 



(2nv) W4(v)dv 



r*+oo 

0.383 k 2 I V 2 C 2 dh 
Jo 



/' A\{f) f- 

00 



' •' X dj\df ; . 
We set x = nfd and find: 

Jr»+co r»+co 
V 2 C 2 d/i J 2 (x) x - s/3 dx 
Jo 



(B.3) 



= 0.360 (V 2 /r Q Y (ro/d) 1/3 = 0.0269 t\ l . 



(B.4) 



The transformation from JB.3I > to ( IB.4I involves a coefficient 
increase by 127r 2 ^ 3 , while the definite integral is equal to 
r(8/3)r(13/6)/[2 8/3 r 2 (ll/6)r(29/6)] = 0.01547. 
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The SF of defocus D^i t) is derived in analogy with the pis- 
ton SF, replacing the response A\( f) for piston with A^if) for 
defocus. The coefficient is 2 times smaller because only one 
aperture is considered. In analogy with jA.8> . 

r*+oo 

D 4 (0 = 0.821 k 2 d 5/i C\ dh K 4 (2tV/d), (B.5) 
Jo 

K A (J3) = 12 [7 3 (x)/x]V 8/3 [l-/oCax)]dje 
Jo 

0.0464 £ 2 + 0.024 /3 6 



1 + 1.2yS 2 +yS 6 



(B.6) 



The approximation has a relative error less than 2% and cor- 
rect asymptotes. Unlike Ki, the Ka function saturates for large 
arguments. Considering only the initial quadratic part of Ka at 
ft <k 1 , we write for small time intervals 

D 4 (f) ~ 0.360 (tV 2 /r ) 2 (r /d) 1/3 = 0.0269 (t/hf. (B.7) 
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